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Abstract 

The complementarity between the quark and lepton mixing matrices is shown to 
provide a robust prediction for the neutrino mixing angle 9^ INS . We obtain this 
prediction by first showing that the matrix Vm, product of the CKM and PMNS 
mixing matrices, may have a zero (1,3) entry which is favored by experimental data. 
Hence models with bimaximal or tribimaximal forms of the correlation matrix Vm 
are quite possible. Any theoretical model with a vanishing (1,3) entry of Vm that is in 
agreement with quark data, solar, and atmospheric mixing angle leads to 6f^ INS = 
(9^2)° • This value is consistent with the present 90% CL experimental upper limit. 
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1 Introduction 



Recent neutrino experiments confirm that the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
[HE] lepton mixing matrix Upmns contains large mixing angles. For example the atmo- 
spheric mixing #23 is compatible with 45° [3] , and the solar mixing 8f 2 MNS is ~ 34° [4] . 
These results should be compared with the third lepton mixing angle 9™ NS which is very 
small and even compatible with zero [HIE], and with the quark mixing angles in the Uckm 
matrix [TJ [8] . 

The disparity that nature indicates between quark and lepton mixing angles has been 
viewed in terms of a 'quark-lepton complementarity' (QLC) [9] which can be expressed in 
the relations 

0PMNS + qCKM „ 45 o . 9 PMNS + qCKM „ (1) 



Possible consequences of QLC have been investigated in the literature [10] and in par- 
ticular a simple correspondence between the Upmns and Uckm matrices has been pro- 
posed QU H21 H3] and analyzed in terms of a correlation matrix [14] [151 [T6| \T7\ [TBI [191 [20] . 
The correlation matrix Vm is simply defined as the product of the CKM and PMNS matri- 
ces, Vm = Uckm • Upmns, and efforts have been done to obtain the most favorite pattern 
for the matrix Vm [20j [21]. Unitarity then implies Upmns — U ckm Vm and one may ask 
where do the large lepton mixings come from? Is this information implicit in the form 
of the Vm matrix? This question has been widely investigated in the literature, but its 
answer is still open (see our section 2). 

Furthermore in some Grand Unification Theories (GUTs) the direct QLC correlation 
between the CKM and the PMNS mixing matrix can be obtained. In this class of 
models, the Vm matrix is determined by the heavy Majorana neutrino mass matrix [T2l [22] . 
Moreover as long as quarks and leptons are inserted in the same representation of the 
underlying gauge group, we need to include in our definition of Vm arbitrary but non 
trivial phases between the quark and lepton matrices. Hence we will generalize the relation 
Vm = Uckm • Upmns to 

Vm = Uckm • ^ • Upmns (2) 

where the quantity Q is a diagonal matrix Q = diag(e 1Wl ) and the three phases u>i are free 
parameters (in the sense that they are not restricted by present experimental evidence). 

The magnitude disparity between the lepton mixing ang le 9f 3 MNS and the other two 
mixings is a rather striking fact. In this paper we carry out the investigation of the 
correlation matrix Vm based on eq. (EI) and prove that it is a zero texture of Vm, namely 
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Vm 13 = 0, that implies a small value for 9f 3 MN with a sharp prediction 

9™ NS = (9±l)°. (3) 

We use the Wolfenstein parameterization for Uckm [23] in its unitary form [24] and param- 
eterize Upmns with the standard phases and mixing angles. As a zero order approximation 
we start inserting by hand the central values of the lepton mixing angles and CKM pa- 
rameters. Owing to the uncertainty in the experimental value for 9± 3 MNS , the possible 
range for the (1,3) entry of matrix Vm may or may not include zero. For example using 
9f 3 M = 3° the (1,3) entry range does not include zero in accordance with eq. (23) in 
ref. [H]. For other choices of 9^ 3 MNS a vanishing (1,3) entry is quite possible, as will be 
seen in section [2j 

It is possible to include bimaximal and tribimaximal forms of the correlation matrix 
Vm in models with renormalization effects [25l |26| [27] that are relevant, however, only in 
particular cases with large tan/3 (> 40) and with quasi degenerate neutrino masses [28] . 
The conclusion for matrix Vm is that the possibility of a bimaximal form, or a tribimaximal 
one is completely open. So in other words, the correlation between the matrices Uckm 
and Upmns is rather nontrivial. 

The investigation we perform for the Vm matrix starts from the fundamental equation 
Vm = Uckm-Q-Upmns and uses the experimental ranges and constraints on lepton mixing 
angles. We resort to a Monte Carlo simulation with two-sided Gaussian distributions 
around the mean values of the observables. The input information on 9i 3 MNS is taken from 
the analysis of ref. [3] which uses neutrino data only. 

The paper is organized as follows: in section [2] we study the numerical ranges of Vm 
entries with the aid of a Monte Carlo simulation, emphasizing on specific points of the 
experimental data. We will show that the vanishing of the (1,3) entry is favored by the 
data analysis. In section [3] we present the matter from a different point of view: we start 
from a zero (1,3) Vm entry (e.g. a bimaximal or tribimaximal matrix) we derive the 
consequent prediction for the Upmns lepton mixing matrix through 

Upmns = (Uckm • ^) 1 ■ Vm (4) 

and the corresponding one for 9f 3 MNS in eq. ([3]). Finally we present a summary and our 
conclusions. 
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2 Which Vm does the phenomenology imply? 



In this section we investigate the order of magnitude of the Vm matrix entries concentrating 
in particular in the (1,3) entry and the important mixing angle 9 13 to which it is directly 
related. We then explicitly study the allowed values of the Vm angles, finally concluding 
that sin 2 6^' = is the value most favored by the data. We will be using the Wolfenstein 
parameterization [23] of the Uckm matrix in its unitary form [21] where one has the relation 

sin 9g KM = A sin 9g KM = AX 2 sin CV FM = AX\p - irj) (5) 

to all orders in A. The lepton mixing matrix Upmns is parameterized in the usual way as 

Upmns = U 23 • $ ■ U 13 ■ $ f ■ U12 ■ $ m . (6) 

Here $ and <3? m are diagonal matrices containing the Dirac and Majorana CP violating 
phases, respectively $ = diag(l, 1, e 1 ^) and $ m = diag(e 1?il , e 1 ^ 2 , 1), so that 



U 



PAINS 



( e^c 12 c 13 e i ^c 13 s 12 s 13 e-^\ 

e l4>1 (-c 23 si2 - e 1 c i2 Si 3 s 23 ) e l<t>2 (c i2 c 23 - e 1 ^ s i2 Si 3 s 23 ) c l3 s 23 
\ e 1?il (— e 1 ^ C12 c 23 si 3 + S12 S23) e 1 * 1 (-e 1<?i c 23 s 12 s 13 - c 12 s 23 ) c 13 c 23 J 



(7) 



2.1 Estimation of Vm entries 

In grand unification models where quarks and leptons belong to the same representation 
of the gauge group, the quark and lepton fields must acquire different phases once their 
symmetry is broken. Hence one should take into account this phase mismatch at low 
energy associated with the form of the CKM and PMNS matrices (0 [7]). To this end we 
introduced the diagonal matrix Q 

Q = dmg{e iuJ1 ,e iLU2 ,e^) (8) 

in the commonly used relational Vm = Uckm • Upmns- This is therefore generalized to 

Vm = Uckm • & • Upmns ■ (9) 

We use for the observed CKM mixing parameters the values A = 0.2237, rj = 0.317, 
p = 0.225, \V cb \ « AX 2 = 0.041, and for the PMNS mixing angles the values 6[ 2 MNS = 34°, 
qpmns _ 450^ qpmns _ 30 pgj For t j ie q phases we resort a Monte Carlo simulation with 



2 see e.g. refs. [9j H4]. 
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flat distributions in the interval [0,27r]. We then get the following range of values for the 
elements of the Vm correlation matrix: 



V, 



M 



I 0.71. ..0.91 0.41. ..0.68 0.10. ..0.22 \ 

0.15. ..0.62 0.40. ..0.74 0.65. ..0.75 
V 0.34.. .0.45 0.54.. .0.64 0.68. ..0.72 / 



(10) 



These values are in good agreement with [14J. The small differences are due to the fact 
that we use the full mixing matrix given in eq. (jHJ) and not the parameterization given in 
eq. (21) of Ref. |14j . Notice that the (1,3) entry of the matrix Vm above cannot be zero, 
so Vm cannot be bimaximal, i.e. of the form 



i 

2 
1 
2 



nor tribimaximal, namely 
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(12) 



where only the absolute values have been considered. The result of eq. (1101) however 
depends on the assumption about the values used for the mixing angles. For example if 
we use a different value for 9[ 3 AINS , namely 6™ NS = 9.2° (see Ref. j3] or our eq. ([TBI for 
the allowed range of 9[ 3 MNS ), we get 



V 



M 



( 0.69. ..0.88 0.39. ..0.67 0.00. ..0.32 \ 

0.09. ..0.67 0.36. ..0.78 0.62. ..0.75 
V 0.28. ..0.51 0.49. ..0.68 0.67.. .0.73 ) 



(13) 



For these values the result is in agreement with the statement that Vm has the (1,3) 
entry equal to zero. It is clear that we need a better investigation of the situation before 
establishing what are the allowed values of the entries of the correlation matrix Vm that 
can be deduced from the experimental data. We next investigate the important entry (1,3) 
as it overwhelmingly affects the 8i 3 MN prediction as will be seen in section [3j 

We parameterize the Vm correlation matrix as the PMNS lepton mixing matrix, i.e. 



Vm = U 23 ■ $ ■ Un -U 



13 



'12 



(14) 



where C/y-s are functions of the mixing angles 0^ M . 
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At first non trivial orders in A we have 



sin 2 # 



2 qV m 
13 



\ _ sin 0PMNS + X gin qPMNS cos qPMNS + Q ^ 



2 



(15) 



It is seen from this expression that the first two terms can cancel each other implying a 
vanishing (1,3) entry of the Vm matrix. In order to better investigate this issue we plot 
in figH] the quantity sin 2 6^' as a function of sin 2 6™ NS . All other observables are fixed 
at their best fit points [3l HI [29] and we allowed the Dirac lepton phase </>, the Majorana 
ones 4>i and 02, and the unphysical phases of Q to vary in the interval [0, 2 7r] with a flat 
distribution. 

As shown in the figure, for the central value of 6™ INS given in [3] the entry (1,3) of 
Vm cannot be zero. However there is a small region (8™ NS « 9.2°) for which 9i£* can be 
zero. This fact has the very important consequence of providing a sharp prediction for the 
unknown mixing angle 6™ N . We will investigate this point in detail in the section [3j 

2.2 The allowed values for tan 2 0^f , tan 2 0$ , and sin^ff 

Here we further investigate the possibility of Vm to be bimaximal or tribimaximal using 
the fundamental equation (Q. We start with a Monte Carlo simulation for the JJ CKM 
parameters, the JJ PMNS mixing angles, the Q and CP phases. 

We use the updated values for the CKM and PMNS mixing matrix, given at 95%CL 
by [29] 

A = 0.2265±g : ggJ2 , A = 0.801183S , 

(16) 



ry = 0.189±°ifi , p = 0.358+ ' 086 



-0.085 j 



with 



\/l - A 2 X 4 (p + if}) 
pj rirj = \L U ■ (17) 



and§ PS] 

sin 2 6™ NS = 0.44 x (l±8:||) , sin 2 9™ NS = 0.314 x (l 
sm 2 9[ 3 MNS = (0.9±g;l) x 10 



+0.18 
0.15 



;is) 



3 The lower uncertainty for sin 2 #13 is purely formal, and correspond to the positivity constraint 
sin 2 6>i3 > 0. 
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With the aid of a MonteCarlo program we generated the values for each variable: for 
the sine square of the lepton mixing angles and for the quark parameters A, A, p, f] we 
took two sided Gaussian distributions with central values and standard deviations taken 
from eqs. (11611181) . For the unknown phases we took flat random distributions in the 
interval [0,2-71"]. We divided each variable range into short bins and counted the number 
of occurences in each bin for all the variables, having run the program 10 6 times. In 
this way the corresponding histogram is smooth and the number of occurences in each 
bin is identified with the probability density at that particular value. A comparatively 
high value of this probability density extending over a wide range in the variable domain 
means a high probability for the variable to lie in this range, therefore that such range 
is 'favoured' by the data being used as MonteCarlo input. Conversely higher probability 
implies better compatibility with experimental data, while lower probability means poor 
or no compatibility with data. 

In figs [2] and [3] we report the results of this simulation. The distributions of tan 2 9^ 
and tan 2 are shown in fig [2j It is seen that the range for which the value of tan 2 9^% is 
compatible with experiments at 90%CL is the interval [0.35, 1.4], so that tan 2 9^ = 1.0 is 
consistent with data. For tan 2 9^ we obtain a range between 0.25 and 1.1 at 90%CL and 
so tan 2 6X2 =1-0 (which corresponds to a bimaximal matrix) only within 3a. Moreover 
the value tan 2 9^ = 0.5 (which corresponds to a tribimaximal matrix), is well inside 
the allowed range. Finally in figE] we plot the distribution for sin 2 9^ . We see that 
sin 2 9^ 3 M = is not only allowed by the experimental data, but also it is the preferred 
value. In the next section we will see that this has important consequences in the model 
building of flavor physics. 



3 Prediction for 0f 3 MNS 

In this section we investigate the consequences of a Vm correlation matrix with zero (1,3) 
entry on the still experimentally undetermined 9f 3 MNS mixing angle. In particular we will 
sec that the 9[ 3 MNS prediction arising from eq. ([9]) or, equivalently, 

Upmns — (Uckm-Q) -Vm (19) 

is quite stable against variations in the form of Vm allowed by the data. 

As previously shown (see section 12.21) . the data favours a vanishing (1,3) entry in Vm- 
So in the whole following analysis we fix sin 2 9^' = 0. We allow the Uckm parameters 
to vary, with a two-sided Gaussian distribution, within the experimental ranges given in 
eq. ffl6l) . while for the Q phases in eq. (jHJ) we take flat distributions in the interval [0, 2ir]. 
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We make Monte Carlo simulations for different values of 9X% and 6^ mixing angles, al- 
lowing tan 2 6X2 and tan 2 9\f to vary respectively within the intervals [0.3, 1.0] and [0.5, 1.4] 
in consistency with the lepton and quark mixing angles (see section 12.21 and fig. [2]). 

In fig JH( left) we plot the distribution of tan 2 8i 2 MNS for values of the correlation matrix 
Vm corresponding to tan 2 9 X2 G {0.3,0.5,1.0} with tan 2 9\ 3 = 1.0. From the figure we 
can check that for tan 2 dXz = 0.3, and 0.5 the resulting distribution for tan 2 6[ 2 MNS is 
compatible with the experimental data. Instead maximal 0X2 and 9\ 3 taken together are 
disfavoured, as the solar angle is hardly compatible with the corresponding allowed interval 
(dot-dashed line). 

In figfjjright) we plot the distribution of tan 2 9 2 ° 3 MNS for tan 2 d\f G {0.5, 1.0, 1.4} with 
tan 0X 2 M = 0.5. Also in these cases we see that the resulting distributions for tan 2 9 23 MNS 
are compatible with the experimental data. 

Finally we report in figj5]the results of our simulation for the quantity sin 2 9{ P 3 MNS . From 
eq. (j!9p . the parameterization of the CKM mixing matrix in eq. (jSJ) and the definition of 
the phase matrix Q in eq. (jSJ) we get 



(U P mns)i3 = e-*" 1 [ ( 1 - ¥.) sin ^e"^ M - A sin 9% cos 9 



13 

+AX\-p + irj + 1) cos 9 v 2 f cos 0[f + 0(A 4 )1 , (20) 



2 f 

iVm „„„ qVm 



so that 



sm 2 9™ NS = sm 2 9Xi'X 2 + 0(X 3 ), (21) 



,2 nPMNS 



where we have used the fact that sin 2 9 13 = and A ps 0(1). 

From eq. f|T9|) and the parameterization used for Vm in eq. (JTlj) we see that sin 2 9f 3 
does not depend on tan 9X 2 M - For this reason the parameter sin 2 9f 3 MN needs to be studied 
as a function of tan 2 9\f only. Fixing for definiteness tan 2 9X2 =0.5 and taking the three 
different values tan 2 9Xf G {0.5,1.0,1.4}, we plot in fig.6 the corresponding distributions 
of sin 2 9f 3 MNS . We note that these values of tan 2 9X 3 practically cover the whole range 
consistent with the data (see fig. 2). 

From fig. [5] it is seen that the sin 2 9[ 3 MNS distributions are quite sharply peaked around 
maxima of 7.3°, 8.9° and 9.8°. Recalling that the shift of this maximum is effectively 
determined by the parameter tan 2 9\ 3 which was chosen to span most of its physically 
allowed range, it is clear that we have a stable prediction for 9± 3 M . 

In order to better clarify this stability, we show in fig. [6] the mean and the standard 
deviation of sin 2 9i 3 MNS obtained with our Monte Carlo simulation for the three chosen 
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values of tan 2 9\ 3 ■ In addition we plot the analytic dependence of sin 2 9[ 3 MNS given by 
eq. (I2ip with the central value of A, the best fit point of sin 2 9[ 3 MNS and its la, 2a and 3a 
from the analysis of ref. [3]. Our prediction for 9i 3 MNS then follows from the experimental 
data on A ,A, p, 77, tan 2 9™ NS and tan 2 9% 3 MNS and the values of tan 2 ^' , tan 2 9^' 
are taken in the intervals [0.3,1.0], [0.5,1.4] respectively, as allowed by the data. For a 
vanishing (1,3) entry of the matrix Vm we finally find 9f 3 MN in the interval [7°, 10°]. 

To conclude this section we note that another prediction for a small 9f 3 MN has recently 
been derived [2QJ 

g PMNS = 9 ° + (sin 3 9g KM ). (22) 

This follows from an assumed bimaximality of a matrix relating Dirac to Majorana neutrino 
states together with the assumption that neutrino mixing is described by the CKM matrix 
at the grand unification scale. Our approach on the other hand is free from any ad hoc 
assumptions. We show that it is a zero texture of the Vm correlation matrix, namely 
Vm 13 = 0, together with all the experimental values of the quark and lepton mixing angles, 
that predicts 9[ 3 MNS = (9±2)°. More importantly we show that the vanishing of this entry 
is favored by the data. Condition Vm 13 = is compatible with Vm being bimaximal (i.e. 
with two angles of 45° and a vanishing one), tribimaximal (i.e. with one angle of 45°, one 
with tan 2 9 = 0.5 and a third vanishing one) or of any other form. Furthermore we make 
use of a phase matrix Q, see eqs. ([319]), that takes account of the mismatch between the 
quark and lepton phases and consider Majorana phases in the Upmns matrix with a flat 
random distribution. 



4 Summary and Conclusions 

In summary, we have investigated the correlation between the CKM quark and PMNS 
lepton mixing matrices, arising in a large class of GUT seesaw models with specific flavor 
symmetries. The detailed analysis developed here uses the fact that the correlation matrix 
is phenomenologically compatible with a tribimaximal pattern, and marginally with a 
bimaximal pattern. This conclusion is different from the one obtained in previous studies 
[T4] and is in agreement with other qualitative arguments that favor the CKM matrix to 
measure the deviation of the PMNS matrix from exact bimaximal mixing |21j . 

In our analysis we found that the mixing parameters tan 2 9X2 and tan 2 9\ 3 vary re- 
spectively within the intervals [0.3,1.0] and [0.5,1.4], while sin 2 9\ 3 l varies in the range 
[0.0, 0.2]. Moreover the preferred value for sin 2 9\ 3 is zero. 

Using these results we investigated the phenomenological consequences of correlation 
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matrices Vm with zero (1,3) entry. The main conclusion of this study is that this large 
class of models is not only compatible with the experimental data, but also that they give 
a robust prediction for 9± 3 MNS mixing angle 



gPMNS = (9 +l )0 (23) 

Whereas the author of ref. [20] obtains a prediction for 9i 3 MNS in a similar range, our result 
cannot be regarded as a straightforward extension or generalization. In fact the condition 
Vm 13 = 0, which is favored by the data, is the only requirement for the prediction ff23l . 
Furthermore we modified the correlation between the CKM and PMNS mixing matrices 
to take account of a phase matrix Q between the quark and lepton fields. Eq. ( |23j) will 
be checked with great accuracy in the next generation of precision neutrino experiments 
(DCHOOZ and others). 

We studied GUT models with flavor symmetry that predict a relation of the type 
Vm = Uqkm ■ • Upmns with Vm 13 = 0. Since in supersymmetric models with tan/3 < 40 
radiative corrections are small [25, 201123 EH], this relation can in such cases be used at low 
energy as in the present paper. Hence if future dedicated experiments exclude 9i 3 MNS ~ 9° 
and supersymmetry is discovered with tan/3 < 40, such models would be ruled out. On 
the other hand, a positive result from 9[ 3 MNS dedicated experiments and tan j3 < 40 would 
be a strong hint for these flavor symmetry models and its specific Higgs pattern. 
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Figure 1: The minimum value allowed for sin 2 OXi' as a function of sin 2 8^ INS . All the 
other CKM and PMNS mixing parameters are fixed at their best fit points given in eq. ffT6l - 
H8l) . The unknown phases u>i, u>2, and u>3 of Q, the Majorana phases fa, and fa, and the 
Dirac one 4> are taken to vary within the interval [0,2 7r] with a flat distribution. We also 
report the values of 9i 3 MNS = 3.0° and 9.2° used in the text. 
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Figure 2: The distributions, normalized to one at the maximum, of tan 2 OX? (solid), and 
tan 2 Q\f (dot-dashed) obtained from the definition of the correlation mixing matrix Vm 
given in eq. (Q by using a Monte Carlo simulation of all the experimental data. 



14 




Figure 3: The distribution, normalized to one at the maximum, of sin 2 0\^' obtained from 
the definition of the correlation mixing matrix Vm given in eq. (jH]) by using a Monte Carlo 
simulation of all the experimental data. We also plot the la and the 2a lines. 




Figure 4: The distribution of tan 2 9[ 2 MNS (left), and tan 2 6f 3 UNS (right) for the CKM 
experimental data and for values of the correlation matrix Vm respectively given by (left) 
tan 2 ^ = 0.3 (dashed), 0.5 (solid), 1.0 (dot-dashed), tan 2 9%$* = 1.0, and sin 2 6^ = 0; 
(right) tan 2 ^f = 0.5 (dashed), 1.0 (solid), 1.4 (dot-dashed), tan 2 9$ = 0.5, sin 2 = 0. 
The shaded areas represent the experimentally allowed regions at 2a for each case. 
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Figure 5: The distribution of sin 2 6i^ INS for the CKM experimental data and for values of 
the correlation matrix Vm given by tan 2 0^ — 0.5, sin 2 Q x ™ = 0, tan 2 Q\f = 0.5 (dashed), 
1.0 (solid), 1.4 (dot-dashed). The shaded area represents the experimentally allowed region 
at la. 
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Figure 6: The allowed values for sin 2 8™ NS as a function of tan 2 O^i 1 under the assumption 
that sin 2 OXz = 0. We report the central and 3<r values obtained from figJHJ and the 
approximate analytical dependence given in eq. (I2TT) . We also plot the experimental central 
value, the lex, the 2<r, and the 3a from [3j. We fixed tan 2 9 X 2 = 0.5 for definiteness. 
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